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Abstract 

Very  large  scale  integration  of  Josephson  junctions  in  a  two-dimensional  series-parallel  ar¬ 
ray  has  been  achieved  by  ion  irradiating  a  YIVhCu  V^  _g  film  through  slits  in  a  nano-fabricated 
mask  created  with  electron  beam  lithography  and  reactive  ion  etching.  The  mask  consisted  of 
15,820  high-aspect  ratio  (20: 1),  35-nm  wide  slits  that  restricted  the  irradiation  in  the  film  below 
to  form  Josephson  junctions.  Characterizing  each  parallel  segment  k,  containing  28  junctions, 
with  a  single  critical  current  Ick  we  found  a  standard  deviation  in  I,±  of  about  16%. 
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In  the  last  two  decades  there  has  been  considerable  effort  aimed  at  developing  a  high-transition 
temperature  ( Tc )  superconductor  Josephson  junction  technology  capable  of  producing  large  num¬ 
bers  of  junctions  with  uniform  electrical  properties,  namely  junction  critical  current  Iq  and  normal 
state  resistance  Rn.1  This  is  especially  challenging  in  high-7j  materials  compared  with  conven¬ 
tional  metallic  low-Tc  superconductors  because  the  superconducting  coherence  length  g  is  much 
shorter  and  highly  anisotropic,  typically  2  nm  in  the  ab  plane  and  0.2  nm  along  the  c-axis  direc¬ 
tion.2  As  a  result,  the  superconducting  order  parameter  is  susceptible  to  structural  and  chemical 
changes  on  atomic  length  scales.  Thus  very  small  imperfections  in  the  Josephson  barrier  or  at  the 
interface  between  the  barrier  and  electrodes  can  drastically  effect  Iq,  since  it  depends  exponentially 
on  the  length  of  the  barrier.  Therefore,  precise  control  at  the  nanometer  scale  is  required  to  make 
multiple  high-rc.  junctions  with  uniform  Iq.  High- 7j  Josephson  devices  are  further  complicated 
by  highly  anisotropic  electrical  transport:  conductivity  along  the  c-axis  is  two  orders  of  magni¬ 
tude  smaller  than  in  the  a-b  plane.3  Such  anisotropy  precludes  the  possibility  of  growing  epitaxial 
multilayers  to  form  sandwich  type  junctions  because  the  highest  quality  thin  films  of  high-Tj.  su¬ 
perconductors  have  c-axes  orientated  normal  to  the  substrate. 

Despite  these  challenges,  a  number  of  junction  fabrication  techniques  have  emerged.  The  earli¬ 
est  was  to  grow  a  thin  film  on  a  bicrystal  substrate  to  create  grain  boundaries  in  the  superconductor 
that  function  as  junction  barriers. 4  These  junctions  are  used  commercially  5,6  in  dc  superconducting 
quantum  interference  devices  (SQUIDs),  which  consist  of  two  junctions  connected  by  a  supercon¬ 
ducting  loop.  While  successful  for  small  scale  applications,  these  junctions  are  not  suitable  for 
large  scale  integration  because  device  layout  is  severely  restricted:  all  the  junctions  must  be  lo¬ 
cated  along  a  single  grain  boundary.  In  addition,  there  are  large  variations  in  Iq  and  Rn  arising  from 
the  short  YBCO  coherence  length  and  imperfections  in  bicrystal  substrates.7  Another  commonly 
used  type  of  junction  is  the  ramp  junction, 8  in  which  the  base  junction  electrode  is  milled  at  an  an¬ 
gle  and  an  intermediate  barrier  material  and  a  superconducting  counter  electrode  are  subsequently 
deposited.  Tens  of  thousands  of  these  junctions  have  been  demonstrated  on  a  single  chip, 9-11  and 
for  some  applications  they  show  considerable  promise  for  a  large  scale  junction  technology. 12,13 
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However,  they  may  not  be  suitable  for  devices  operating  at  liquid  nitrogen  temperatures,  because 
the  operating  temperature  is  typically  below  50  K. 14  Furthermore,  the  inter-junction  spacing  can 
only  be  scaled  down  to  a  few  gm, 15  which  does  not  allow  them  to  be  used  in  high  frequency 
devices  that  may  require  junction  spacings  of  the  order  of  100  nm. 16 

An  alternative  technology  is  to  use  ion  damage  Josephson  junctions,  which  are  fabricated  by 
using  ion  irradiation  to  reduce  Tc 17  in  a  narrow  region  of  a  superconducting  bridge.  Above  the 
narrow  region  Tc  and  below  the  undamaged  electrode  Tc,  the  irradiated  area  behaves  as  a  Joseph¬ 
son  junction.  Since  ion  damage  junctions  have  no  interfaces  between  different  materials,  they 
can  be  arbitrarily  positioned,  and  densely  spaced  (~  100  nm). 18  One  way  of  forming  this  type  of 
junction  is  to  pattern  a  thick  mask  to  protect  the  superconducting  electrodes  and  to  damage  un¬ 
masked  regions  with  ion  irradiation.  This  technique  was  first  demonstrated  by  Tinchev  to  fabricate 
SQUIDs. 19  Others  have  followed  with  variations  of  the  process  to  reproduce  single  junctions,20-24 
series  arrays  of  tens  of  junctions,25,26  and  a  series  array  of  280  SQUIDs.27  Here  we  report  the 
scaling  up  of  this  process  by  two  orders  of  magnitude,  with  a  2-dimensional  (2D)  array  of  15,820 
(28  x  565)  Josephson  junctions  or  equivalently  15,255  (27  x  565)  SQUIDs. 

We  designed  a  SQUID  configuration  with  incommensurate  loop  areas  because  of  recent  sig¬ 
nificant  interest  in  these  structures.  The  critical  currents  of  each  SQUID  oscillate  with  incommen¬ 
surate  periodicities  as  a  function  of  applied  magnetic  field  B.  As  a  result,  the  critical  current  of 
the  entire  array  Ica  attains  a  maximum  value  at  zero  magnetic  field,  but  is  nearly  zero  elsewhere. 
Consequently,  for  Josephson  junctions  with  nonhysteretic  current- voltage  (I-V)  characteristics,  the 
array  biased  with  a  (negative)  current  exhibits  a  sharp  peak  in  voltage  at  zero  field.  This  behav¬ 
ior  was  first  demonstrated  by  Sohn  et  al.2&  using  2D  arrays  of  Nb-A10rNb  Josephson  junctions. 
Later  Carelli  et  al 29  showed  that  one  could  use  a  series  array  of  incommensurate  area  SQUIDs 
as  an  absolute  magnetometer.  More  recently,  arrays  of  incommensurate  area  SQUIDs  made  from 
the  high-7).  superconductor  YBa2Cu307_g  (YBCO)  have  been  fabricated  with  loops  connected  in 
series,  parallel,  and  a  series-parallel  combination. 30,31  It  has  been  suggested  that  it  may  be  possible 
to  use  incommensurate  area  SQUID  arrays  as  radio  frequency  (RF)  amplifiers. 32 
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The  layout  of  our  array  is  shown  in  Figure  1 .  For  efficient  coupling  of  RF  in  future  experiments 
we  chose  a  microstrip  line  configuration  with  SQUID  loops  cut  into  it.  The  following  criteria  were 
used  in  designing  this  array.  The  width  of  the  microstrip  was  set  to  be  500  jUm  so  that  its  impedance 
would  be  50  Q.  on  a  0.5-mm  thick  sapphire  substrate.  The  junction  widths  were  chosen  to  be  2  /tm 
to  avoid  self-field-limited  junction  effects33  and  to  keep  Rn  high.  The  parameter  /3/  =  2LIq / <t>o 
(where  L  is  the  loop  inductance  and  <t>o  is  the  flux  quantum)  was  kept  small  ( /3/  <  0.05)  to  ensure 
that  Ica(B)  modulates  almost  to  zero.34  This  was  achieved  by  using  loop  areas  <  90  /im2  to  keep 
L  <  53  pH.  The  overall  length  of  the  array  in  the  direction  of  the  bias  current  was  chosen  to  be 
3  mm,  less  than  1/4  wavelength  for  frequencies  below  10  GHz,  so  that  the  array  behaves  as  a 
lumped  element.  The  500-/im  width  of  the  microstrip  enabled  us  to  fit  27  loops,  connected  in 
parallel,  across  it  (Figure  1).  The  widths  of  these  loops  ranged  from  9.5  to  22.5  /im  in  0.5-/im 
increments.  In  the  direction  of  the  current,  565  parallel  segments  with  lengths  ranging  from  3.0  to 
4.0  /im  in  0.25-/im  increments  were  connected  in  series  in  113  blocks  of  five. 

To  create  this  array,  a  200-nm  thick  YBCO  thin  film35  was  thermally  coevaporated  on  a  sap¬ 
phire  wafer  followed  by  a  gold  contact  layer  deposited  in  situ.  We  patterned  the  films  using  pho¬ 
tolithography  and  Ar+  ion  milling  to  fabricate  the  microstrip  with  15,255  SQUID  loops  etched  into 
it  (Figure  2a).  The  gold  layer  over  the  junctions  was  removed  using  a  subsequent  photolithography 
step  and  Transene  chemical  gold  etch,  leaving  the  contact  pads  (shown  in  Figure  1).  The  next 
step  was  to  fabricate  the  junctions.  The  wafer  was  coated  with  a  700-nm  layer  of  Shipley  SI 808 
hard-baked  photoresist  that  served  as  the  main  ion  stopping  layer.  A  25 -nm  layer  of  germanium 
was  electron-beam  evaporated  on  top  of  the  resist  and  served  as  an  etch  stop.  We  next  spun  100 
nm  of  polymethylmethacrylate  (PMMA)  resist  on  to  the  Ge  for  electron-beam  lithographic  pat¬ 
terning  (Figure  2b).  Using  a  100-keV  Leica  VB6-HR  nanowriter,  we  wrote  35-nm  wide  lines  in 
the  PMMA  over  the  locations  intended  for  the  junctions.  This  pattern  was  transferred  into  the  Ge 
layer  with  reactive  ion  etching  (RIE)  in  a  HBr-Cl2  plasma.  The  pattern  in  the  Ge  was  transferred 
to  the  resist  using  low  temperature  (-100°C),  low  pressure  (5  mtorr)  oxygen  RTF  (Figure  2c).  This 
process  resulted  in  a  high  aspect  ratio  of  the  line  width  to  the  trench  depth  (1:20)  because  resist 
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etches  two  orders  of  magnitude  faster  then  Ge  in  an  oxygen  plasma. 

Figure  3a  shows  a  scanning  electron  microscope  (SEM)  image  of  a  section  of  the  structure 
after  etching  and  before  implantation.  The  image  was  taken  using  an  18-kV  acceleration  voltage 
that  partially  penetrated  the  mask  and  allowed  for  simultaneous  imaging  of  both  the  top  of  the 
mask  and  the  YBCO  film  below.  A  magnified  view  of  two  junction  regions  (Figure  3b)  shows  the 
5-jUm  long,  35-nm  wide  slits  patterned  into  the  mask.  The  etches  were  checked  for  undercutting 
by  simultaneously  processing  a  silicon  dummy  wafer  with  the  same  layered  structure.  This  test 
wafer  was  cleaved,  and  the  cross-sectional  image  is  shown  in  Figure  3c.  Following  etching,  the 
wafers  were  implanted  commercially  with  200-keV  Ne+  at  a  dose  of  lxlO13  ions/cm2.  The  Tc  of 
the  YBCO  not  protected  by  the  mask  was  lowered  by  ion  damage,  forming  planar  inline  Josephson 
junctions. 

Completed  devices  were  attached  to  a  printed  circuit  board  equipped  with  ^-filters  and  a  silicon 
diode  thermometer.  The  board  was  mounted  inside  a  vacuum  probe  and  cooled  in  a  liquid  nitrogen 
bath.  We  measured  the  resistance  of  the  device  as  a  function  of  decreasing  temperature  using  a 
lock-in  amplifier  and  a  22-Hz,  8-/1A  peak  current.  As  shown  in  Figure  4,  the  device  exhibited 
two  distinct  superconducting  transitions:  one  for  the  undamaged  YBCO  and  another  for  the  ion 
damaged  material  (junctions).  The  transition  temperatures  were  86.3  ±0.2  K  and  78.7  ±0.2  K, 
where  we  defined  each  transition  temperature  to  correspond  to  the  point  of  inflection  in  the  R  vs. 
T  characteristic,  or  equivalently  to  the  peaks  in  dR/dT  vs.  T .  Thus,  the  ion  damage  reduced  Tc 
by  about  7.6  K.  The  full-width  at  half  maximum  (FWHM)  for  the  dR/dT  peak  of  the  irradiated 
material,  1.4  ±0.2  K,  was  similar  to  that  of  the  electrodes,  1.3  ±0.2  K.  Thus,  the  spread  in  the  Tc 
of  the  junctions  in  the  565  segments  in  series  is  close  to  the  spread  in  the  Tc  of  the  unirradiated 
YBCO.  This  implies  that  the  spread  in  Tc  along  each  parallel  segment  is  likely  to  be  comparable. 
These  measurements  suggest  that  the  ion  damage  is  quite  uniform  across  the  array,  implying  that 
the  resistances  of  the  weak  links  are  similarly  uniform.  We  note  that  we  previously  observed  a 
similar  spread  in  the  Tc  of  a  ID  array  of  280  SQUIDS.27 

We  measured  I-V  characteristics  for  the  array;  data  for  five  representative  temperatures  are 
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shown  in  Figure  5.  To  fit  the  characteristics,  we  made  three  initial  assumptions.  First,  we  assumed 
that  each  parallel  segment  has  a  resistance  of  R/565,  where  R(T)  is  the  fitted  resistance  of  the  total 
array  at  each  temperature.  Second  we  assumed  that  each  parallel  segment  k  of  the  array  behaves  as 
a  single,  resistively  shunted  junction36,37  with  a  critical  current  Ick.  Third,  we  assumed  that  the  set 
of  Ick  have  a  normal  distribution  given  by  Ick  =  Ic(\  +  Sk) ,  where  Ic  is  the  mean  critical  current  of  the 
parallel  segments  and  8k  is  the  fractional  variation  in  the  critical  current  Ick  of  segment  k.  Under 
these  assumptions,  we  fitted  the  data  using  a  nonlinear  regression  with  V  =  (R/565)Y,k=i(I2  ~ 
l2k)  '/2  (7  >  Ick)  to  find  Ic,  and  the  standard  deviation  in  Ic,  o(Ic).  We  found  that  cr(/c)  converges 
to  about  16%  with  decreasing  temperature.  At  higher  temperatures  thermal  noise  rounding38  of 
the  I-V  characteristic  was  indistinguishable  from  spreads  in  Ic.  In  the  inset  of  Figure  5,  we  see  that 
the  resistance  determined  from  the  fit  increases  rapidly  with  increasing  temperature.  The  solid  line 
linear  fit  shows  that  Ic  scales  as  (1  —  T / Tc )3,  consistent  with  a  soft  boundary  model  of  ion  damage 
junctions  near  Tc.  We  remark  that  in  this  2D  array  the  parallel  segments  enable  us  to  operate 
with  a  lower  critical  current  in  each  junction  and  retain  phase  coherence  closer  to  Tc  where  the 
soft  boundary  model39  is  applicable.  In  our  previous  work  on  a  one-dimensional  series  array27  the 
need  to  achieve  phase  coherence  required  us  to  work  further  below  Tc  where  the  critical  current 
scales  as  (1  —  T /Tr)2.  We  found  that  r(/c)=78.3  K  for  Tc  =  8  /iA,  in  good  agreement  with  the 
transition  temperature  found  from  Figure  4  at  the  bias  current  of  8  /iA  used  for  that  measurement. 

To  study  the  response  of  the  array  to  a  magnetic  field,  the  array  was  biased  with  a  static  (neg¬ 
ative)  current  and  the  voltage  across  the  array  was  measured  as  a  function  of  applied  magnetic 
field.  Three  different  field  scales  are  shown  in  Figure  6.  For  large  fields  (Figure  6a)  the  Fraunhofer 
diffraction  pattern  (approximately  |sinc  (/r<t>/<t>o)  |)  from  the  areas  of  the  individual  junctions  is 
observed.  Five  oscillations  are  visible,  indicating  that  the  junction  areas  have  good  uniformity. 
From  the  period  in  the  Fraunhofer  pattern  we  calculate  the  YBCO  penetration  depth  to  be  480 
nm  at  77.3  K.  The  hysteresis  in  the  data  arises  from  vortex  trapping  that  depends  on  the  magnetic 
field  sweep  direction  and  magnitude.40  For  intermediate  fields  (Figure  6b)  both  the  Fraunhofer  and 
SQUID  pattern  are  visible.  Detailed  measurements  near  B  =  0,  for  6  different  bias  currents,  are 
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shown  in  Figure  6c.  The  central  peak  is  observed  at  38  pT.  This  is  where  our  applied  field  cancels 
the  component  of  Earth’s  field  perpendicular  to  the  plane  of  the  array  leaving  it  in  zero  absolute 
field.  Away  from  the  peak,  the  oscillations  diminish  as  a  result  of  destructive  interference  due  to 
the  different  loop  sizes.  We  performed  simulations  to  characterize  the  width  of  the  central  voltage 
peak,  A B,  which  we  define  as  the  separation  of  the  minima  on  either  side.  From  Figure  6c,  we  find 
A B  =  20.5  ±  0.5  T.  We  inserted  the  areas  of  the  SQUIDs  and  the  fitted  parameters  I,  and  R  into 
Eq.  (14)  of  Oppenlander  et  al ,41  The  results  of  the  simulations  greatly  underestimated  the  width 
of  the  peak.  We  found,  however,  that  a  modified  model,  in  which  the  27  x  5  block  of  SQUIDs  was 
modeled  by  a  ID  series  array  of  5  SQUIDs  each  with  areas  corresponding  to  the  average  areas  of 
the  parallel  segments  (48,  52,  56,  60,  64  ^m2),  predicted  A B  -  20.9  T,  about  2%  greater  than  the 
measured  value.  We  have  no  ready  explanation  for  this  empirical  result. 

We  also  observe  a  linear  tilt  in  V  vs.  B,  visible  in  Figure  6c,  that  we  ascribe  to  the  motion  of 
Josephson  vortices  in  the  array.  This  motion  is  analogous  to  the  motion  of  Abrikosov  vortices  in 
type  II  superconductors.  Here,  the  bias  current  exerts  a  Forentz  force  given  by  (7  x  7?)  where 
~7  is  the  current  density  and  ~B  is  the  magnetic  field.  This  force  causes  vortices  to  drift  with  a 
steady  velocity  across  the  array,  creating  an  electric  field  in  the  direction  of  the  bias  current.  On 
our  measurements  this  appeared  as  magnetoresistance.  For  bias  currents  in  the  opposite  direction, 
the  asymmetry  had  the  opposite  sign,  consistent  with  our  interpretation.  We  analyzed  our  data 
using  the  Bardeen-Stephen  equation42  p  =  ^>qB/t]B,  which  relates  the  flow  resistivity  p  to  the 
magnetic  field  B,  to  determine  J],  the  vortex  viscosity  coefficient.  Using  the  bias  current  and 
junction  geometrical  parameters  we  converted  V  to  p,  and  determined  d>o B/rj  from  a  linear  fit  of 
the  data.  We  estimate  p  =  3  x  1(U7  kg  m_1s_1  at  a  current  bias  of  100  pA.  We  are  not  aware 
of  any  other  measurements  of  Josephson  vortex  motion  in  2D  arrays  of  this  type  for  comparison. 
However  this  value  is  typical  for  the  viscosity  coefficient  of  Abrikosov  vortices  in  high-transition 
temperature  superconductors. 

In  conclusion,  by  building  an  array  of  15,820  Josephson  junctions,  we  have  demonstrated  that 
it  is  possible  to  fabricate  large  numbers  of  ion  damage  high-7j  Josephson  junctions  with  uniform 
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properties  over  large  areas.  This  technique  shows  promise  for  a  very  large  scale  junction  tech¬ 
nology.  In  particular,  the  relatively  small  standard  deviation  (16%)  in  I(±  implies  that  long,  linear 
arrays  containing  parallel  segments  of  equal-area  SQUIDs  have  substantial  potential  for  applica¬ 
tions  such  as  magnetometers  and  amplifiers.  Since  the  flux-to-voltage  transfer  coefficient  scales  as 
N ,  the  number  of  SQUIDS,43  and  the  signal-to-noise-ratio  scales  as  A1/2,  such  arrays  are  likely 
to  overcome  the  low  transfer  coefficient  and  signal-to-noise  ratio  associated  with  the  low  junction 
resistance  inherent  in  this  technology.  Furthermore,  smaller  barrier  dimensions  combined  with  nar¬ 
rower  bridges  made  with  refinements  in  lithography  and/or  new  techniques  involving  nano-wires44 
may  allow  for  substantial  increases  in  IqRu ,  opening  up  this  process  to  more  applications  such  as 
rapid  single  flux  quantum  logic  or  precision  digital-to-analog  conversion. 
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Figure  1:  Layout  of  a  2D  array  of  SQUIDs  patterned  as  a  microstrip.  The  array  has  565  series 
segments  of  28  Josephson  junctions  connected  in  parallel.  Junctions  are  represented  by  an  x. 
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Figure  2:  Fabrication  process  for  ion  damage  Josephson  junctions,  (a)  Magnified  view  of  a  pat¬ 
terned  YBCO  bridge,  (b)  Layered  structure  of  photoresist,  germanium  and  PMMA.  (c)  Structure 
after  EBL,  RIE  and  during  irradiation,  (d)  YBCO  with  reduced  Tt  shown  in  red. 
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Figure  3:  SEM  photographs  after  RIE  etching  of  the  implant  mask,  (a)  A  segment  of  the  device 
showing  some  of  the  YBCO  loops,  (b)  Magnified  view  over  two  of  the  slits  that  were  used  to 
pattern  the  junctions,  (c)  Cross-sectional  view  of  a  dummy  sample  to  test  the  RIE  etch  for  under¬ 
cutting  and  to  ensure  that  the  etching  was  complete. 
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Figure  4:  Resistance  R  and  dR/dT  vs.  temperature,  (a)  Transitions  of  the  electrodes  and  weak 
links,  (b)  Transition  of  the  weak  links  shown  on  an  expanded  scale.  The  width  of  the  peak  in 
dR/dT  at  each  transition  gives  a  measure  of  the  homogeneity  of  Tc  of  the  electrodes  and  weak 
links. 


Figure  5:  TV  characteristics  for  5  temperatures  for  the  2D  Josephson  array  showing  both  data 

-  1  /3 

(circles)  and  fitted  results  (solid  lines).  The  inset  shows  fitted  R  vs.  '/’(red)  and  Ic  vs.  '/"(blue). 
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Figure  6:  Voltage  as  a  function  of  magnetic  field  for  (a)  large  ,  (b)  medium,  and  (c)  small  field 
ranges.  The  smallest  range  shows  measurements  for  six  bias  currents. 
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